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Abstract 



Diamond and zinc-blende photonic crystals are studied both in the purely 
dielectric case and in the presence of small inclusions of a low absorbing 
metal. It is shown that small metal inclusions can have a dramatic effect on 
the photonic band structure. Several complete photonic band gaps (CPBG's) 
can open in the spectrum, between the 2nd-3rd, 5th-6th, and 8th-9th bands. 
Unlike in the purely dielectric case, in the presence of small inclusions of 
a low absorbing metal the largest CPBG for a moderate dielectric constant 
(e < 10) turns out to be the 2nd-3rd CPBG. The 2nd-3rd CPBG is the 
most important CPBG, because it is the most stable against disorder. For a 
diamond and zinc-blende structure of nonoverlapping dielectric and metallo- 
dielectric spheres, a CPBG begins to decrease with an increasing dielectric 
contrast roughly at the point where another CPBG starts to open-a kind 
of gap competition. A CPBG can even shrink to zero when the dielectric 
contrast increases further. Metal inclusions have the biggest effect for the 
dielectric constant e € [2,12], which is a typical dielectric constant at near 
infrared and in the visible for many materials, including semiconductors and 
polymers. It is shown that one can create a sizeable and robust 2nd-3rd 
CPBG at near infrared and visible wavelengths even for a photonic crystal 
which is composed of more than 97% low refractive index materials (n < 1.45, 
i.e., that of silica glass or a polymer). In the case of silica spheres with 
a silver core, the 2nd-3rd CPBG opens for a metal-volume fraction f m ~ 
1.1% and has a gap width to midgap frequency ratio of 5% for f m 2.5%. 
Within the 2nd-3rd CPBG of 5%, absorption remains very small (< 2.6% 
once the CPBG is centered at a wavelength A > 750 nm), which should 
be tolerable in most practical applications. The metallo-dielectric structures 
display a scalinglike behavior, which makes it possible to scale the CPBG from 
microwaves down to the ultraviolet wavelengths. Aluminum, copper, and gold 
cores yield almost identical results, provided that sphere radius r s > 250 nm. 
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For r s < 250 nm the results for different metals can be increasingly different 
with decreasing r s , nevertheless, qualitative features remain the same. These 
findings open the door for any semiconductor and polymer material to be 
used as genuine building blocks for the creation of photonic crystals with 
a CPBG and significantly increase the possibilities for experimentalists to 
realize a sizeable and robust CPBG in the near infrared and in the visible. One 
possibility is a construction method using optical tweezers, which is analyzed 
here. 

PACS numbers: 42.70.Qs, 78.67.-n, 82.70.Dd, 87.80.Cc 
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I. INTRODUCTION 



Photonic crystals are structures with a periodically modulated dielectric constant. In 
analogy to the case of an electron moving in a periodic potential, light propagating in a 
photonic crystal experiences multiple scattering leading to the formation of Bloch waves and 
photonic band gapsBHi. If the band gap persists for both polarizations and all directions of 
propagation one speaks of a complete photonic band gap (CPBG)Bi. Light with frequencies 
within a CPBG is totally reflected since it cannot propagate inside the crystal. Yet light can 
propagate through waveguides carefully designed within such a photonic crystal, even when 
such a waveguide has sharp bends!. In the last decade, photonic crystals have enjoyed a 
lot of interest in connection with their possibilities to guide light and to become a platform 
for the fabrication of photonic integrated circuitsi'i. There is a common belief that, in the 
near future, photonic crystals systems will allow us to perform many functions with light 
that ordinary crystals do with electrons. In addition to numerous potential technological 
applications (filters, optical switches, superprisms, cavities, etc!), photonic crystals also 
promise to become a laboratory for testing fundamental processes involving interactions of 
radiation with matter under novel conditions^!. The presence of a CPBG causes dramatic 
changes in the local density of states, which offers the possibility to control and engineer the 
spontaneous emission of embedded atoms and molecules!. 

Unfortunately, technological difficulties in fabricating CPBG structures rapidly increase 
with decreasing wavelength for which a CPBG is required. Despite the research activities 
of a large number of experimental groups, achievement of a CPBG below infrared (IR) 
wavelengths for both two- and three-dimensional (3D) photonic structures is still elusive, 
mainly because the required dielectric contrast 5 to open a CPBG is rather high!. Here 
S = ^max/^min, where e max (smin) is the dielectric constant of a material component, used to 
fabricate a photonic crystal, with the largest (smallest) value of e. For the best geometries 
5 ~ 5 is required!!!. Already this threshold value of 5 excludes the majority of semiconduc- 
tors and other materials, such as (conducting) polymers, from many useful photonic crystal 
applications. However, for applications one needs a sufficiently large CPBG to leave a mar- 
gin for gap-edge distortions due to omnipresent defects!. Let us define g w as the gap width 
to midgap frequency ratio, Aw/w c . Then in order to achieve g w larger than 5%, 5 > 9.8 
and 5 ^ 12 is required for a diamond! and face-centered-cubic (fee) structure!, respectively. 
This leaves only a couple of materials for photonic crystals applications at near infrared and 
optical wavelengths^!. 

Surprisingly enough, it will be demonstrated here that there is a way to create a sizeable 
and robust CPBG in the near infrared and in the visible even for a photonic crystal which 
contains more than 97% of a material with the refractive index below n = 1.45, i.e., that 
of silica glass or a polymer. The trick is to place small inclusions of a low absorbing metal 
in the right dielectric structure. It turns out that the effect of small inclusions of a low 
absorbing metal on the photonic band-gap structure is very strong for a parent diamond 
dielectric structure. In the latter case, the actual metal-volume fraction f m needed to open 
a CPBG of more than 5% depends on the available material dielectric constant e and can be 
kept below 1% for e s = 4 (Fig. | below). Given a metal- volume fraction f m , the strongest 
effect on CPBG was found when both spheres in the lattice primitive cell contained identical 
metal cores. Surprisingly, the inclusions have the biggest effect for the dielectric constant 
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e G [2, 12], which is a typical dielectric constant at near infrared and in the visible for many 
semiconductors and polymers. Given the desired gap width to midgap frequency ratio of 5%, 
the smallest absorption was found when only one of the two spheres in the lattice primitive 
cell contained a metal core, i.e., for a zinc-blende structure, even though the required f m 
was typically twice as large as that for a diamond structure. This result suggests that, 
given the metal filling fraction f m , absorption is reduced due to a reduction of near-field 
electromagnetic energy transfer between the metal cores with an increased separation of the 
metal islands in the structure. For a zinc-blende structure, absorption within a CPBG of 
5% can be kept below 2.6% once the CPBG is centered at a wavelength A > 750 nm. The 
results on the absorption are by far the best which have been demonstrated for a 3D metallo- 
dielectric structure with a CPBG. They are an order of magnitude better than for the case 
of an fee lattice of metal coated spheres0. Our results on the absorption compare well to 
the best results for one- dimensional and two-dimensional metallo-dielectric photonic crystals 
which show absorptance of ~ 1% and ~ 3%, respectively, at A ~ 600 mxiS. Photonic band- 
structure calculations also revealed a surprising scalinglike behavior of our metallo-dielectric 
diamond and zinc-blende structures, which is only intrinsic to purely dielectric structures. 
The scalinglike behavior means that once a CPBG is found, the CPBG can be open for any 
wavelength, simply by scaling all the sizes of a structure-an extremely useful property from 
a practical point of view. 

A metal-core dielectric-shell sphere morphology seems to play an essential role in the 
effect of small metal inclusions on photonic band structure. If the same volume of a metal is 
spread homogeneously within the spheres, no CPBG is found. Similarly, a diamond structure 
of small metal nanospheres embedded in a dielectric matrix shows a much smaller CPBG 
and much higher absorptance for a comparable f m than either a diamond or a zinc-blende 
close-packed structure of metal-core dielectric-shell spheres in air. On purely experimental 
grounds, only the case of spheres with a metal core is investigated here. Usually, a metal 
shell around a dielectric core is formed by an aggregation of small metallic nanoparticles. 
The shell has to be around 20 nm thick before it becomes completeEl. With an emphasis on 
photonic structures in the near infrared and in the visible, the 20 nm shell thickness then 
would mean a rather high threshold value of the metal filling fraction f m (of the order of 
5%). On the other hand, it is much easier to tune the metal filling fraction f m from zero to 
a few percent by coating small metal nanoparticles with a dielectric in a controlled way0. 
Moreover, a dielectric shell is necessary to prevent aggregation of the metallic particles by 
reducing the Van der Waals forces between them. In the latter case, a coating of roughly 
20 nm is required. Having metal cores also means that an active light-emitting dielectric 
material (semiconductor or polymer) is not isolated as in the case of metal shells. Hence 
crystals of metal@dielectric particles can perform more easily some useful optoelectronic 
functions. 

Before proceeding with the photonic band-structure calculation of metallo-dielectric di- 
amond and zinc-blende structure of spheres, it was necessary to recalculate the photonic 
band structure of the diamond lattice of dielectric spheres. The previous plane wave method 
(PWM) calculations of the photonic band structure of diamond^ and zinc-blended lattices 
of dielectric spheres were found to be plagued by large errors. Therefore, after describing 
the method of our calculation in Sec. [H], we discuss and analyze the differences with respect 
to earlier photonic band-structure calculations for a diamond lattice of dielectric spheres by 
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Ho, Chan, and Soukoulisi. The effect of small inclusions of a low absorbing metal on the 
photonic band structure of a diamond and a zinc-blende structure is investigated in Sec. [II. 
In Sec. [1 11 A| the dependence of CPBG's on the metal volume fraction f m is investigated. 
The dependence of gaps on either the shell dielectric constant in the diamond case, or on 
the dielectric constant of the dielectric sphere in the zinc-blende case, is summarized in Sec. 



I [ IB} . A surprising scalinglike behavior of the metallo-dielectric structures is then discussed 
in Sec. [Ill CJ Reflectance, transmittance, and absorptance of light incident on the metallo- 
dielectric structures is dealt with in Sec. I IIDj . Most of the results will be shown for the 



case of silver cores. However, as discussed in Sec. [IV A| , where aluminum, copper, and gold 
cores are examined, the results so obtained are almost insensitive to the type of metal used, 
provided that the sphere radius is larger than « 250 nm. If the sphere radius is smaller, 
most of the results are still qualitatively valid. We then end up with a discussion in Sec. [TV] 
and conclusions in Sec. [V|. 



II. BENCHMARKING DIAMOND STRUCTURE OF DIELECTRIC SPHERES 

Indisputably, one of the hallmarks in the study of photonic band gap structures is an 
article by Ho, Chan, and Soukoulisi. It opened a way for the fabrication of the first photonic 
structure with a complete photonic band gap (CPBG)i and advanced the field considerably. 
Its main conclusion is that, regarding a CPBG, a diamond structure of (overlapping) spheres 
fares much better than a simple fee structure: (i) a CPBG opens between the 2nd and 3rd 
bands (the 8th- 9th bands for an inverted fee structure), and, consequently, is much more 
stable against disorder!, (ii) the threshold value of the dielectric contrast 5 to open the 
CPBG is 4 (8.2 for an inverted fee structured), and (iii) the CPBG is significantly larger 
(15% and 5%, for the respective diamond and fee close-packed lattices of spheres with a 
dielectric contrast 5 = 12.96). 

Results of our calculations are summarized in Figs. Jl], |2|. They were obtained using the 
photonic Korringa-Kohn-Rostocker (KKR) method§0'0"(see also Appendix A). The KKR 
method can be used for scatterers of arbitrary shapefflj and is optimized for lattices of 
spheres. The photonic band gap structure was calculated with the value of the angular- 
momentum cutoff Imax = 9. Stability of the Ewald summation was thoroughly checked. 
Convergence of the method with increasing l max is demonstrated in Fig. |^. The two lowest 
bands in Fig. [I] (but not the higher ones) agree well with those calculated in Ref. 0. In 
agreement with Refs.BEl, no CPBG was found for the case of a diamond lattice of nonover- 
lapping air spheres in a dielectric. However, for the inverse case of a diamond lattice of 
nonoverlapping dielectric spheres in air, the earlier results of Ho, Chan, and SoukoulisS on 
the photonic band structure were found to be nonconverged. Take, for example, the test 
case of a close-packed diamond lattice of nonoverlapping dielectric spheres with the dielectric 
constant e s = 12.96 in air (Fig. [I]). A nonconvergence of the results of§ is already witnessed 
by an unusually large deviation of the effective refractive index n e ff, as calculated from the 
band structure in the L direction, from n^ff = y/s e ff, where e e ff is calculated by the Garnett 
formula!!, 

e eff ne h (1+2 fa)/(l -fa), (1) 
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which usually provides a very good fit to n e //. Here / is the sphere volume fraction (/ = 
0.34 for closed-packed diamond and our zinc-blende lattices of spheres) and a is a sphere 
polarization factor [a = (e s — Sh)/{e 8 + 2e/i), where e s and Eh are the dielectric constant 
of sphere and of host medium, respectively]. In our case, the effective refractive index 
n e ff = 1.604 is only 10% larger than n^f = 1.456. For comparison, the upper edge of the 
bands at the L point in Ref.i seems to be almost 20% lower than in our case, resulting 
in more than a 30% deviation from 'nMf. Not surprisingly, 6 years after Ref.@ the PWM 
(see Fig. 4 of Refill) yielded for the same test case a 2nd-3rd CPBG width of only « 8%, 
half of that presented in RefR In the meantime, Soziier, Haus, and Ing uvai, p. 13 971, 
again using the PWM, argued that the 2nd-3rd CPBG is less than 3.5%, and they called 
for a recalculation of their results using more efficient methods, such as the photonic KKR 
method. 

For a diamond lattice of nonoverlapping dielectric spheres, i.e., for sphere filling fraction 
f s varying from to the close-packed case f cp = 0.34, two CPBG's can occur simultaneously 
(see Fig. p]), between the 2nd-3rd bands, and, as in an inverted fee case, between the 
8th-9th bands!. According to Fig. |2|, the photonic band structure of a diamond lattice of 
nonoverlapping dielectric spheres shows a much more complex behavior of CPBG's than an 
fee lattice. For an fee lattice of air spheres in a dielectric, the calculated gap width to midgap 
frequency ratio g w of the 8th-9th CPBG is monotonically increasing with the increasing 
dielectric contrast Eh/£ s and shows a saturated behavior for different filling fractions (see 
Fig. 4 in Ref.0). A remarkable feature of a diamond lattice of nonoverlapping dielectric 
spheres is a kind of gap competition. According to Fig. 0, a CPBG first increases with the 
dielectric contrast e s /eh, but, at some dielectric contrast shortly before another CPBG starts 
to open, the CPBG begins to decrease and it may eventually disappear. As the dielectric 
contrast e s /eh increases further, the first CPBG may reappear again, with its reappearance 
being signaled by a decrease of the second CPBG. Fig. |^ illustrates that kind of behavior 
for f a = 0.17 and f s = f cp . Contrary to previous calculations^, the 2nd-3rd CPBG is found 
not to be the dominant one. It satisfies the bound of 3.5% by Soziier, Haus, and Inguva 
(see Ref.0, p. 13 971) not only for e s ~ 13 but apparently for all e s . For f s = f cp , the 
2nd-3rd CPBG only persists for e s G [5.2, 16.3] and does not exceed 2.3% (for e s — 9). The 
2nd-3rd CPBG is then absent till e s = 40, where it reappears again. We disagree with Ref.0 
that the 2nd-3rd CPBG vanishes for / < f cp . It is true that, for e s = 12.96, the 2nd-3rd 
CPBG closes already for f s = 32%, i.e., for f s only 2% smaller than the close-packed case. 
However, with increasing e s the 2nd-3rd CPBG can appear at smaller and smaller filling 
fractions. For instance, for e s ~ 20 the 2nd-3rd CPBG begins to open for f s = 17%, but 
again its g w does not exceed 2.9%, satisfying the 3.5% boundEil. The dominant CPBG is the 
8th-9th CPBG. For e s = 12.96, the 8th-9th CPBG persists down to f s = 4%. For /, = 17% 
and e s = 12.96 it can reach 12%, and it reaches its maximum of ~ 14% for e s = 15.5. For 
fs — fcp it reaches a maximum of ~ 14.7% for e s = 35. However, the threshold value of e s 
for its opening is 7.9, comparable to that for an inverted fee latticeS. Hence, regarding a 
CPBG, a diamond lattice of nonoverlapping dielectric spheres fares still better than an fee 
structure of air spheres in a dielectric^, but the advantage is much smaller than has been 
previously thought. 

An argument to persuade the reader about the correctness of our result is that they can 
be reproduced, up to some minor differences, by the PWM based MIT ah initio programil. 
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However, unlike the case of a simple lattice (one scatterer per lattice primitive cell)l3, for 
the case of a diamond lattice of dielectric spheres, even when using the PWM based MIT 
ab initio program!!, one has to take a much larger number of plane waves than expected to 
reach a convergence comparable with the photonic KKR method. As shown by Megensil, to 
reach convergence of the photonic band structure of a diamond lattice of dielectric spheres 
within 1%, the number of plane waves N c has to exceed 32 768 (cf. Ref.S) and still an 
extrapolation iV — > oo0 has to be performed. Hence, the iV — > oo extrapolation0 might 
also be necessary to ensure a convergence for the MIT ab initio programed when modeling 
systems with sharp material discontinuities for other complex lattices. Note that N c is two 
orders of magnitude larger than Nhr* = 339 which reproduces the results of§ for the whole 
range of filling fractions < / < It turns out that some other results computed in the 
early days of the plane-wave expansion method can be imprecise and plagued by large errors 
also, an example being Ref.il. 



III. EFFECT OF METAL INCLUSIONS ON PHOTONIC BAND STRUCTURE 

In this section, the effect of small inclusions of a low absorbing metal on the photonic 
band-gap structure for a parent close-packed diamond and zinc-blende dielectric structure of 
spheres is investigated. The close-packed case of spheres in air was chosen because this is one 
of the most probable cases to be fabricated either by "do-it-yourself organization" il or by a 
microrobotic technique^. At the same time, even when starting from colloidal suspension^, 
one expects to use dried structures in air for practical applications which will most probably 
consist of touching spheres (close-packed structures). A metal core will be placed inside 
either both dielectric spheres in the lattice primitive cell, or only one single sphere in the 
lattice primitive cell. In the former case, the resulting structure will be a diamond structure, 
whereas in the latter case the structure becomes a zinc-blende structure. In the following 
we only consider the special case of zinc-blende structures of spheres, when both spheres in 
the lattice primitive cell have identical radii. 

To reach convergence of bands with 1% precision in the metallo-dielectric case, the pho- 
tonic KKR methodQOEa with the angular momentum cutoff parameter l max = 8 was used. 
The KKR method is best suited to deal with highly dispersive scatterers. Computational 
time with and without the dielectric constant dispersion is the same and only depends on the 
angular momentum cutoff parameter (see also Appendix A). As in earlier work0!!!!, only 
the real part of the material dielectric constant was used when the photonic band structure 
was calculated. Reflectance, transmittance, and absorptance were calculated by adapting 
available computer co which is based on the layer photonic KKR (LKKR) methodic. 
In the latter case, both the real and imaginary part of the material dielectric constant were 
used. The dielectric constant of a metal was chosen according to Ref.@. Above the vacuum 
wavelength of ~ 1900 nm, Palik's data for aluminum, copper, and gold are no longer avail- 
able. In this case a Drudelike fit to a metal dielectric constant was used with parameters 
determined by Ordal et a/.!! (Ordal's parameters has also been used by El-Kady et a/.0 to 
fit a metal dielectric constant from far infrared down to optical wavelengths.) The same 
procedure for the band structure and absorption calculations has also been employed in 
an earlier work on metallo-dielectric structures^!!!!. It is true that a complex dielectric 
constant turns Bloch eigenvalues into complex resonances in the lower-half of the complex 
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frequency planeE-i However, it can be shown that the imaginary part of a low absorbing 
metal has onlylittle influence on the projections of the complex resonances on the real 
frequency axisSZI. This can also be shown by a direct comparison of the "projected" band 
structure against reflection and transmission calculations (see Sec. p 1 [T3|) . Note in passing 
that a comparison of theorvand experiment using the photonic KKR methods has so far 
been more than satisfactory!!'!!. 

Photonic band-structure calculations demonstrated that several CPBG's can open in 
the spectrum. Depending on the metal filling fraction f m < 10% and a typical dielectric 
constant e G (1,20] of a dielectric at near infrared and in the visible, a diamond structure 
can have CPBG's between the 2nd-3rd and 8th-9th bands, whereas a zinc-blende structure 
can have CPBG's between the 2nd-3rd, 5th-6th, and 8th-9th bands!!. Figure f| shows the 
photonic band structure for a close-packed zinc-blende lattice of spheres in air. One of the 
two spheres in the primitive lattice cell is a silver core - n s = 1.45 (silica) shell sphere with 
r c/ r s — 0.75 (f m = 7.2%), whereas the other is a homogeneous dielectric sphere with £ — 12 
of the same radius r s = 80 nm. The frequency is in units [c/A], as in Fig. [l|. By comparing 
Fig. [I] and Fig. || it is obvious that as the result of the metal inclusions much larger 2nd-3rd 
and 8th- 9th CPBG's open. There is also an additional CPBG between the 5th and 6th 
bands. Metal inclusions seem to align band edges at the L and W points and the degenerate 
band edge of the 3rd and 4th bands at the L point raises above the band edge of the 3rd 
band at the U point. 

In the following subsection, the dependence of the gaps on the metal-volume fraction is 
investigated. The dependence of gaps on either the shell dielectric constant in the diamond 
case, or on the dielectric constant of the dielectric sphere in the zinc-blende case, is summa- 



rized in Sec. [j 1 1B[ A surprising scalinglike behavior of the metallo-dielectric structures is 
then discussed in Sec. [Ill C| . Reflectance, transmittance, and absorptance of light incident 
on the metallo-dielectric structures is dealt with in Sec. [IIID| . 



A. Dependence of gaps on the metal-volume fraction 

Photonic band-structure calculations revealed a strong effect of metal inclusions on the 
CPBG between the 2nd and 3rd bands. Given a metal-volume fraction f m , the largest 2nd- 
3rd CPBG opens for a close-packed diamond lattice of metal core-dielectric shell spheres in 
air (Fig. |[), followed by a close-packed zinc-blende structure (Figs. ||, |7]). Unlike in the 
purely dielectric case, which has been discussed in Sec. OL in the presence of small inclusions 
of a low absorbing metal the dominant CPBG for a diamond lattice is that between the 2nd 
and 3rd bands (Figs. |Sp. The 2nd-3rd CPBG is the most important CPBG, because it is 
the most stable against disorder!. In the zinc-blende case, the 2nd-3rd CPBG is also much 



larger than in the purely dielectric case and is a dominant CPBG for e G (1,9] (Figs. 10, 
il~D . For higher e the dominant CPBG becomes the 5th-6th CPBG (Figs. [H], [II]). Photonic 
band-structure calculations revealed several remarkable features of the diamond and zinc- 
blende metallo-dielectric structures. First, a strong increase of the 2nd-3rd CPBG with f m , 
once a threshold value of f m is reached. Figure § shows that if the silver cores are placed 
inside spheres with a refractive index n s = 1.45 (silica) and a radius of 80 nm, the 2nd-3rd 
CPBG below 600 nm opens for f m « 1.1% and reaches 5% already for f m « 2.5%. If the 
sphere refractive index n s increases further and approaches the threshold refractive index 
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contrast of ~ 2.3, for which the 2nd-3rd CPBG of the parent diamond structure of non- 
overlapping dielectric spheres begins to open (see Fig. H), the required metal / m 's rapidly 
decrease. For n s = 2 (ZnS) spheres with a silver core, the 2nd-3rd CPBG begins to open 
already for r c /r s = 0.13 (f m ~ 0.07%) and reaches 5% already for r c /r s ps 0.26 (f m ~ 0.6%). 
Figure ^| shows that in the zinc-blende case one needs almost twice as large f m as in the 
diamond case to obtain a comparable 2nd-3rd CPBG. For example, to obtain the 2nd-3rd 
CPBG of 5% in the case when silver cores are only placed inside a single n s = 1.45 (silica) 
sphere in the lattice primitive cell, one needs f m ~ 5%. Figure |7| shows the effect of the 
dielectric sphere morphology on the 2nd-3rd CPBG, which is either a n c = 2 (ZnS) core - 
silica shell dielectric sphere with fixed r c /r s = 0.60 (dashed line) or a homogeneous silica 
sphere (solid line). The other sphere in the lattice primitive cell is a silver core - n s = 1.45 
(silica) shell sphere. 



B. Dependence of gaps on the dielectric constant 

Figures [5], [5], |T0 , and [TTJ illustrate the effect of a varying dielectric constant on a CPBG. 
In the diamond case, the shell dielectric constant is varied. In the zinc-blende case, one of 
the spheres in the lattice primitive cell remains a silver core - n s = 1.45 (silica) shell sphere 
with fixed r c /r s = 0.75 (f m = 7.2%), and we vary the dielectric constant of the dielectric 
sphere. Quite amazingly, small metal inclusions have the biggest effect for e G [2, 12] (cf. 
Figs. [10] and |11]), which is a typical dielectric constant in the near infrared and in the 
visible for many materials such as semiconductors and polymers. The dependence of the 
CPBG's on the dielectric constant is very complex. Figures |^, [lTj, and [TTJ show that the 



2nd-3rd CPBG remains the largest CPBG for e < 9. According to Figs. [10] and [TTJ, for 
zinc-blende structures and e > 9, the 2nd-3rd CPBG becomes smaller than the 5th-6th 
CPBG. The 8th-9th CPBG remains small for both diamond and zinc-blende structures. As 
already seen for a diamond structure of non-overlapping dielectric spheres (see Fig. 0)1, 
if e increases beyond a certain threshold value, the 2nd-3rd CPBG begins to narrow and 
may even disappear. Such a narrowing of the 2nd-3rd CPBG with an increasing dielectric 
constant seems a characteristic feature of the diamond and zinc-blende structures. Figure 



TTJ shows the same kind of "gap competition" which has been first observed in the purely 
dielectric case (see Fig. ^]): a CPBG first increases with increasing dielectric constant, but, 
at some dielectric contrast shortly before another CPBG starts to open, the CPBG begins 
to decrease and it may eventually disappear. 



C. Scalinglike behavior 

A further remarkable feature of the metallo-dielectric structure is a surprising scalinglike 
behavior, which is intrinsic only to ideal dispersionless structures. Given the metal filling 
fraction f m , the best scalinglike behavior is observed in our zinc-blende structures (see Figs. 



T2| , |T3| , |T4l , |T7] ). A scalinglike behavior is expected once the metal e at a given midgap 
wavelength becomes sufficiently large and negative. Indeed, on the single-scatterer level it 
makes rather little difference if the real part of the metal e = —200 or e = — oo, the limit of 



a perfect metal. Figure 15 shows that differences in the sphere extinction efficiency rapidly 
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decrease with an increasing sphere radius. There is only a little difference in the extinction 
efficiencies of r s = 300 nm and r s = 600 nm spheres. One can, therefore, view the preceding 
figures involving r s = 300 spheres as a limiting case of a perfect metal. For silver inclusions, 
the midgap position and the 2nd-3rd CPBG edges follow almost perfect scaling down to 
r s = 200 nm. For r s < 200 nm, the gap width and gap edges begin to depend stronger 
and stronger on decreasing r s . However, if the midgap wavelength is plotted in nm against 
sphere radius, one would see a linear dependence down to r s = 80 nm (see Figs. O, ITBI). 
for which the 2nd-3rd CPBG midgap wavelength is typically below 600 nm and the real 
part of silver e ~ —10. Therefore, in Fig. O the ratios of the midgap and gap edges of 



the 2nd-3rd CPBG to the sphere radius are plotted to make deviations from the scalinglike 
behavior visible at all. This scalinglike property is very useful from a practical point of view. 
It means that once a CPBG is found, with some midgap wavelength A c , the CPBG can be 
centered at any other wavelength by a simple scaling of all the sizes of a structure. 



D. Reflectance, transmittance, absorptance 

In order to further investigate optical properties of our photonic structures, reflectance, 
transmittance, and absorptance were calculated. By reflectance, 71, and transmittance, 
T, we mean total reflected and transmitted field (the specular and all propagating higher 
diffraction orders), respectively. Absorptance, A, is then defined as the total loss in the 
structure, A = 1 — 1Z — T. Reflectance, transmittance, and absorptance were calculated by 
adapting available computer co which is based on the photonic LKKR methodic. The 
same method was also used in Refs.0 andS To reach convergence within 1% around the 
2nd-3rd CPBG, l max = 6 was used. 

For a metallo-dielectric structure is characteristic that, beyond a certain crystal thickness, 
reflectance, transmittance, and absorptance within a CPBG change only negligibly with an 
increasing crystal thicknessS As a function of the crystal thickness, £, absorptance behaves 
as 

A(£) = A sat [i-ex P (-e/e A )}, (2) 

where £a is an "absorption" length and A sat is the saturated absorptance. For some CPBG's, 
A sat can only be a few per cent. This saturated behavior can easily be explained as follows. 
Because of a CPBG, light can only penetrate the first few layers of a crystal. Hence, an 
addition of further layers beyond the penetration length has only minor influence. As a direct 
consequence of the saturation, absorptance per crystal thickness decreases approximately 
linearly with the crystal thickness, A{£)/£ ~ A sa t/£- For a high metal- volume fraction, 
as for an fee photonic crystal, such a saturated behavior can already be established after 
three crystal planes0H. In our case, the metal volume fraction is almost 20 times smaller. 
Consequently, the saturation (with two digit precision) is established after 12 crystal planes 
for a normal incidence in the (111) crystal direction, or after two unit-cell thickness. The 
(111) crystal direction has been chosen because this choice makes it possible to compare the 
properties of our photonic crystals against those of fee photonic crystals, which properties 
have been studied in the same crystal directionQH. It is useful to remind that the (111) 
direction corresponds to the T-L direction in the notation of special points on the surface 
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of the Brillouin zone. Hence, for our diamond and zinc-blende structures, reflection and 
transmission measurements in the (111) direction probe the slice of band structure (see 
Figs. [I] and ^) which consists of the line parallel to the y axis which intersects the x axis at 
the L point. 

In Fig. [T7| reflectance, transmittance, and absorptance of light incident normally on a two 
unit cells (12 planes) thick zinc-blende lattice of spheres in air stacked in the (111) direction 
is shown. Note first a scalinglike behavior of reflectance, transmittance, and absorptance 



in Fig. O. The spheres are as in Figs. 12, O, [L4[ The AABBCC stacking sequence is 



depicted in Fig. ITS]. The gap width to midgap frequency ratio of the 2nd-3rd CPBG of 



these configurations has been shown in Fig. O. The midgap wavelength of the 2nd-3rd 



CPBG has been shown in Figs. [13], [14|. 

Calculations based on the LKKR methodffl showed that the saturated absorptance 
within a CPBG can be kept at very small levels. Given the desired gap width of 5%, 
the smallest absorption was found for close-packed zinc blende structures, even though the 
required f m was typically twice as large as that of the close-packed diamond structure. 
According to Fig. |T7|, the saturated absorptance within a CPBG of > 5% can be kept below 
5% for a CPBG centered at A ~ 600 nm and is less than 2.6% for a CPBG centered at 
A > 750 nm. This should be tolerable in most practical applications. The results on the 
absorption are by far the best which have been demonstrated for a 3D metallo-dielectric 
structure with a CPBG. They are an order of magnitude better than for the case of an fee 
lattice of metal coated spheres0. (Although we have used the same method as in Ref.B§ to 
calculate absorption, when recalculating the results in Ref.0 for an fee lattice of metal-coated 
dielectric spheres, a much larger absorptance was obtained.) Our results for the saturated 
absorptance compare well to the best results for one-dimensional (ID) and two-dimensional 
(2D) metallo-dielectric photonic crystals which show a saturated absorptance of ~ 1% and 
~ 3%, respectively, within a CPBG centered at A ~ 600 nm@. 

Note in passing that within the 2nd-3rd CPBG and for normal incidence in the (111) 
direction, almost 99% reflection goes into the specular beam. As the angle of incidence 
increases, the specular reflectivity decreases and gradually more light is reflected in non- 
specular directions. 



IV. DISCUSSION 

A. Use of others metals as a sphere core 

So far, when presenting our results on the effect of small inclusions of a low absorbing 
metal on the photonic band structure of diamond and zinc-blende photonic crystals we have 
used exclusively silver. This noble metal guarantees the best_properties of metallo-dielectric 
photonic crystals in the near infrared and in the visi bleSiH. However, almost identical 
results have also been obtained for aluminum, copper, and gold cores, provided that sphere 



radius r s > 250 nm (see Fig. |19|). For r s < 250 nm spheres, the results for different metal 



can be increasingly different with decreasing r s , nevertheless, qualitative features remain the 



same (see Fig. |20|) . 



This behavior is explained as follows. When the scalinglike behavior was discussed in 



Sec. [HI C[ it was argued that one can view the results involving r s > 250 nm spheres as those 
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corresponding to a limiting case of a perfect metal, when metal dielectric constant is e = — oo. 
Indeed, Fig. 15 showed that differences in the single-sphere extinction efficiency rapidly 



decrease with an increasing sphere radius. Therefore, one expects only minor differences for 
diamond and zinc-blende photonic crystals with r s > 250 nm spheres when gold, copper, or 
some other metal is used instead of silver, with the other lattice and sphere parameters being 
identical. These expectations are fully confirmed by exact calculations, and the results are 
displayed in Figs. 0, 120, A similar behavior, i.e., only a slight change in absorption when 
gold or copper has been used in place of silver, has also been found for an fee structured. 



B. Fee vs diamond and zinc-blende structures 

Obviously, small metal inclusions do not open a CPBG in every dielectric structure. For 
example, a simple fee lattice of dielectric spheres with a metal core requires f m ~ 50% to 
open a CPBGc§B3. For an fee lattice of metal-coated dielectric spheres, the required metal 
filling fraction f m is slightly lower but still very high (« 4O%0) and still almost 40 times 
larger than in our case. Therefore, not surprisingly, when going further to shorter and shorter 
wavelengths, one is facing an increasing absorption: at A ~ 600 nm the absorption exceeds 
30% even within a CPBG. Although such a metallo-dielectric fee structure could provide 
a at near infrared, the extension to the visible is difficult. A probable reason 

of such a big effect of small metallic inclusions on the photonic band gaps of the diamond 
and zinc-blende structures compared to an fee structure is because the former have much 
better photonic band-gap properties already in the limiting cases of purely dielectric and 
purely metallic spheres. In the limit of dielectric spheres in airiB, the threshold value of the 
dielectric contrast to open the 2nd-3rd CPBG for a diamond structure is only 5.2@, whereas 
an fee structure does not have any CPBG, irrespective of the sphere dielectric const antfE. 
In the opposite limit of pure metallic spheres, the 2nd-3rd CPBG for the diamond lattice is 
huge. For silver, depending on the sphere radius r s , it can stretch from 60% (r s = 80 nm) 
to 75% (r s > 300 nm)iil. This is consistent with a previous estimate of g w > 60% for the 
case of an ideal metal (e s = — oo)0. On the other hand, for an fee lattice a higher CPBG 
opens (between the 5th-6th bands) and is smaller 40%)E3. 



C. Behavior of the 2nd-3rd CPBG with an increasing metal-volume fraction 

In order to explain the opening of a CPBG with an increasing metal volume fraction, the 
following argument has been used by Zhang et aiS Low-frequency waves propagate through 
a 3D metallo-dielectric structure as in an effective medium characterized by the effective 
dielectric constant £ e //& This low-frequency pass band will have zero group velocity when 
the wave vector approaches the Brillouin-zone boundary where the wave vector k « 7r/a, a 
being the lattice constant. Thus the highest frequency oji of this low-frequency pass band 
~ c/ (ay/£ e ff). On the other hand, metal cores inscribe a void with a scale of the same 
order as a, and the fundamental resonance modes in the void should have frequencies of 
the order uo u Pd c/a. These voids are always connected together by channels, so that the 
modes can hop from one void to another forming a pass band. Since e e ff increases when the 
metal-volume fraction increases, it may be possible to have uoi < uo u when the metal-volume 
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fraction exceeds a certain threshold value. Then, assuming that there are no other bands 
between the low-frequency pass and hopping bands, a CPBG should be formed between ui 
and uj u . However, Zhang et a/.H applied their argument to an fee lattice where a CPBG 
opens between the 5th-6th bands, i.e., the band on the lower edge of the CPBG is not a 
low frequency pass band. It seems more natural to reconsider their argument for our case 
of diamond and zinc-blende photonic crystals of coated spheres. 

The Garnett formula@ [Eq. ([]])] yields e e // for a cubic lattice of both homogeneous and 
coated spheres. One has only to substitute a corresponding polarization factor for a. In the 
case of a coated sphere with a single coating, embedded in the host medium characterized by 
the dielectric constant e^, let us denote the dielectric constant and radius of the sphere core 
(shell) by e c and r c (e a and r s ), respectively. Defining x = r^/r^, a c = (e c — e s )(e c + 2e s ), and 
a s = (e s — Sh){s s + the polarization factor a of a coated sphere with a single coating 
can be written a si 

a s + xa c (e h + 2e s )/ (e s + 2e h ) 

a = ■ (3) 

1 + 2xa c a s 

Here < x < 1 and < ctj < 1. Note that a c (a s ) coincides with the polarization factor 
of a homogeneous sphere with the dielectric constant e c (e s ) embedded in the host medium 
characterized by the dielectric constant e s (Polarization factors of a homogeneous 

sphere are also recovered from a in one of the following limits: x — > 0, x — > 1, and e c — > £ s @). 

For a diamond lattice of close-packed metal-core dielectric-shell spheres in air, e e ff, as 
calculated from Eqs. ([]]) and @, increases from « 1.3 up to ~ 1.9, i.e., by only a factor of 
« 1.5, when r c /r s increases from 0.1 to 0.75 (irrespective of the real part of the dielectric 
constant of the metal core when varied between —50 and —500, since the metal-volume 
fraction is kept rather small). Therefore, in our case, the change in s e ff with an increasing 
metal-volume fraction f m is by no means dramatic. Yet the behavior of the 2nd-3rd CPBG 
edges, as shown in Fig. |2T], seems to be qualitatively in agreement with the argument by 
Zhang et a/H For r s = 80 nm spheres, both CPBG edges first decrease with increasing 
r c /r s . However, when r c /r s increases beyond ~ 0.5, the upper 2nd-3rd CPBG edge slows 
down its decrease and becomes approximately constant, whereas the lower edge of the 2nd- 
3rd CPBG continues in its monotonic decrease. The case of r s = 300 nm spheres shows a 
qualitatively different behavior of the upper 2nd-3rd CPBG edge, which frequency increases 
with increasing r c /r s . This behavior is exactly what one would expect according to the 
argument by Zhang et a/.c2l. In the case of r s = 300 nm spheres, the 2nd-3rd CPBG is 
located at the frequency range for which silver skin depth [= c/(2u Im n{uj))\ is only slightly 
above 10 nm, whereas the silver core has radius between 150-180 nm. Therefore, light 
can only penetrate a small fraction of the silver core. Consequently, since the void volume 
between spheres shrinks as r c /r s increases, the fundamental resonance modes in the void 
shift to higher frequencies. In the case of r s = 80 nm spheres, the 2nd-3rd CPBG is located 
at the frequency range for which silver skin depth can be as large as 50 nm, i.e., light can 
penetrate entire silver core. Therefore, the "void" part of the argument by Zhang et al. is 
much less applicable in the r s = 80 nm case. 
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D. Coated versus homogeneous metal spheres 



We have also examined if a large and robust CPBG can open by using homogeneous 
metal nanospheres arranged on a diamond lattice in a dielectric, with the possibility of 
fabricating a diamond structure by layer-by-layer deposition. According to Fig. [Tj|, the 
sphere projections on the stacking z direction in neighboring planes do not overlap if the 
sphere radius r s < r cp /3, where r cp is the sphere radius in the close-packed case. Therefore, if 
r s < r c P /3 it would be possible to fabricate a diamond lattice by depositing first a hexagonal 
plane of spheres on a substrate and then filling in the interstitial by a dielectric untill all 
spheres are embedded in a dielectric layer. One would then repeat the procedure layer by 
layer and grow gradually a diamond lattice. Unfortunately, the constraint r s < r cp /3 means 
that the resulting sphere filling fraction f s < 1.26%, which is too small to achieve a large 
and robust CPBG with homogeneous metal nanospheres in a dielectric. Indeed, no CPBG 
was found for a diamond structure of homogeneous metal nanospheres with r s < 0.45r cp 
(fm ~ 3.1%) embedded in a dielectric host with the dielectric constant e — 2.1, 3, 4, and 
6.25. Photonic band-structure calculations revealed that one needs f s = f m ~ 10% to open 
the 2nd-3rd CPBG of 5% with metal nanospheres in silica. For comparison, note that for 
a close-packed diamond structure of n s = 1.45 dielectric spheres with a metal core and 
radius of 80 nm, the 2nd-3rd CPBG opens for r c = 0.32r cp , and for r c = 0.45r cp the 2nd-3rd 
CPBG reaches almost 5% (Fig. ^). This suggests that a metal core-dielectric shell sphere 
morphology plays a crucial role in the formation of CPBG's. Indeed, if the same volume 
of a metal is spread homogeneously within the spheres, no CPBG opens in the spectrum. 
In the latter case, the sphere dielectric constant was calculated using the Garnett formula 

A dielectric contrast between the dielectric shell and the host dielectric medium also 
seems to be an important factor. Figure [22] shows the dependence of the 2nd-3rd CPBG for 
a close-packed metallo-dielectric diamond lattice of spheres of radii r s = 80 nm and 300 nm 
on the host dielectric constant. Spheres have a silver core and dielectric shell with e = 8, 
with the core to total sphere radii ratio either r c /r s = 0.5 (f m = 4.25%), or r c /r s = 0.6 
(fm — 7.3%). The shell dielectric constant has to be sufficiently large to open a CPBG. For 
example, for the dielectric shell with e = 4 and r c /r s = 0.4 (f m = 2.2%), the 2nd-3rd CPBG 
is closed in the range of the host dielectric constant in Fig. 1221. 



E. Absorption 

Similarly, as in the case of a photonic crystal of small r s = 5 nm Drude plasma spheres^, 
one observes (see Figs. [T7|, [19], |20|) that the absorptance peaks coincide with the transmit- 
tance peaks. An explanation for this behavior can already be traced down to the case of 
a homogeneous slab of metal. In the latter case, one can show that neither the imaginary 
part of the metal dielectric constant nor the inverse skin depth (absorption coefficient) char- 
acterize properly the absorption of an incident light. Indeed, the absorption coefficient is a 
measure of absorption of only that light which has entered the metal. Higher transmittance 
means higher probability of light being absorbed. Hence, a proper measure of the degree of 
absorption of an incident light is rather a product of the transmittance and the absorption 
coefficient. 
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In the purely dielectric caseEj and in the case of an fee crystal of small r s = 5 nm Drude 
plasma spheres with sufficiently large relaxation time of conduction electrons and a low 
volume fraction (/ = 0.1)0, transmittance peaks occur when the effective half wavelength, 
A/2 = Tc/k±, k± being the wave vector in the stacking direction, times an integer equals 
approximately the thickness of the photonic crystal slab. Hence, when transmittance is 
plotted against k±, one observes a series of almost equidistantly spaced peaksSS. It is well 
known that dk±/du can attain very high values near band edges. As a rule of thumb, band 
edges are located either at the surface, or in the center (k = or T point) of the Brillouin 
zone. At the points on the surface of the Brillouin zone with reflection symmetry, frequency 
eigenvalues satisfy oo(k) = cj(— k). For those points dk^/du even diverges to infinity. (The 
latter is equivalent to saying that at such a band edge the group velocity of light approaches 
zero.) Therefore, if transmittance is plotted against frequency, one observes a considerable 
increase in the density of transmittance peaks in the immediate vicinity of a frequency 
gap. However, when the relaxation time of conduction electrons of Drude plasma spheres 
gradually decreases (what is equivalent to an increasing imaginary part of the dielectric 
constant), the densely spaced sharp transmittance and absorptance peaks near band edges 
merge and generate a rather smooth profiled. This is what is also observed in our case, as 
exemplified in Figs. [T7|, and |2D|, which show no densely spaced and sharp transmission 
peaks in the proximity of the lower 2nd-3rd CPBG edge. 

A dramatic change in the absorptance at the lower and upper banclgap edges is a 
reminiscent of the Borrmann effect, previously observed in x-ray scatterings. The origin of 
this effect is that fields in the proximity of band-gap edges approach a standing wave. In the 
proximity of one band-gap edge, fields are mostly localized in the regions of a low dielectric 
constant, whereas in the proximity of the other band-gap edge, fields are mostly localized 
in the regions of a high dielectric constant. This sort of behavior can already be observed 
for a simple ID periodic Bragg stack of dielectric layers^. Figures [H], [19], and ^0| suggest 
that in the proximity of the lower band gap edge of our zinc-blende structures electric field 
is mostly localized in air, whereas in the proximity of the upper band gap edge electric field 
is mostly localized at sphere positions, leading to a significantly higher absorption. 

Our results on absorption compare well to the best results for ID and 2D structures^!. 
In these respective cases, the absorption of an optimized metallo-dielectric structure can 
be ~ 1% and ~ 3% within a CPBG centered at A ~ 600 nm. Given the metal filling 
fraction f m , it has been found that absorption is reduced when metal forms isolated islandsEx 
This is an expected result: when metal forms an interconnected network, DC conductivity 
is nonzero and long-range conduction currents are induced which lead to higher losses. 
Our results suggest that absorption is reduced further when separation of metal islands 
increases. Indeed, a larger separation of metal cores in our zinc-blende structure compared 
to a diamond structure resulted in a significantly reduced absorption. This is probably 
because of a reduction of near-field electromagnetic energy transfer between the metal cores 
with an increased separation of the metal islands in the structure. 



F. Fabrication 

Whereas a simple fee structure occurs naturally in colloidal crystals formed by 
monodisperse colloidal suspensions of microspheres or via template-directed colloidal 
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crystallization!!, diamond and zinc-blende structures are very difficult to fabricate. Nev- 
ertheless, several new ways have recently been proposed to make these and other complex 
structure and therefore diamond and zinc-blende structures deserve our full attention. 

One of the possibilities is a "do-it-yourself organization" which, using artificial templates, 
enables one to create structures beyond thermodynamic limits^. Patterning of a surface 
can be performed using optical tweezersEl and particle manipulation can also be employed 
in a microrobotic technique^. Difficulties in optical trapping of a metal-core dielectric-shell 
sphere increase with an increasing core size. There is always a danger that a metal core can 
be melted in the laser focus. For spheres with a 7.5 nm gold core and a silica shell with 
a final radius of 79 nm, optical trapping has been recently demonstrated experimentally in 
RefH In Fig. |23| we show that metal-core silica-shell spheres with much larger cores can 
also be optically trapped. We plot there the axial trapping efficiency Q for silver core - 
n s = 1.45 (silica) shell spheres with fixed r c /r s = 0.3 immersed in water (n^ = 1.33). The 
axial trapping efficiency Q, or the dimensionless normalized axial force, when multiplied by 
the focused laser beam power P, determines the resulting axial force F on a sphere immersed 
in a homogeneous medium with refractive index nh, 

F = (n h /c)PQ. (4) 

The laser wavelength A is increased by 100 nm with increasing sphere radius, from A = 600 
nm for r s = 80 nm to A = 1000 nm for r s = 200 nm. Q is plotted against q/r s , the center 
offset from the focus in units of the sphere radius. The beam opening angle of the focused 
beam was taken to be 6 = 78 and the ratio of the objective focal length to the beam waist 
was set to 0.99. The normalized axial force Q was calculated by adapting results of Ref.El. A 
stable equilibrium location on the beam axis corresponds to the point where the normalized 
axial force Q, considered as a function of q/r s , crosses zero with a negative derivative. Fig. 
ZB"| shows that the spheres can be optically trapped for all radii, with a stable equilibrium 



at around q/r s pa 2. Note that the equilibrium point is almost five times further from the 
focus in the current case than in the case of purely dielectric spheres (cf. Fig. 2 of Ref.0). 
A larger distance from the focus is convenient experimentally as it reduces the chance of 
melting a metal core by an intense light. 



G. Related work 

In our study we have only considered metal-core dielectric-shell particles which can easily 
be fabricated^. In this respect, our results regarding a CPBG are not optimized. Probably 
the required metal-volume fraction can be rendered even smaller than in our case. One way 
is to use perturbation theory along the lines of Ref.@. Here it has been demonstrated that 
a photonic band gap can be enlarged by placing metal insertions at particular positions of 
a latticed. However, this work dealt with an enlargement of already existing band gaps, 
whereas, in our case, small inclusions of a low absorbing metal are used to open a CPBG. 
Also Ref.@ only deals with photonic crystals in two dimensions, whereas our work covers 
3D photonic crystals. Photonic band gaps of AB(3) and B-3 structures of metallo-dielectric 
spheres have been investigated in Ref.c3. It has been shown that an AB(3) photonic crystal 
contains a large photonic band gap and a B-3 photonic crystal has more than one gap when 
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the filling ratio of metal spheres exceeds a threshold. A complementary case of 3D diamond 
structures of metal-coated dielectric spheres has recently been discussed in Reffi Such 
metal nanoshells, consisting of a dielectric core with a metallic shell of nanometer thickness 
(cca 50 nm), can be used to vary the optical resonances of such nanoparticles over hundreds 
of nanometers in wavelength by varying the relative dimensions of the core and shelfi It 
has been shown that, under certain conditions, it is possible to center a resonance band, 
which originates from the single-sphere Mie resonances of metal-coated dielectric spheres, 
at the middle of the 2nd-3rd CPBG. The position and width of the resonance-induced in- 
gap modes can be controlled by making appropriate choices of the dielectric constant of the 
inner dielectric sphere and thickness of the metal-coating layer, respectively. Such properties 
can be very useful in making optical band filters as well as microcavity lasers if they are 
sustained in the presence of dissipation. In the case of metal-core dielectric-shell spheres 
the tunability of the single-sphere Mie resonances is much smaller. Nevertheless, in the 
latter case the resonance-induced in-gap modes can be observed for certain diamond and 
zinc-blende configurations. A study of these resonance-induced in-gap modes in the case of 
metal-core dielectric-shell spheres will be presented elsewhere. Note that in order to center 
a single-scatterer resonance band within a CPBG, the presence of a metal is crucial. The 
latter is impossible in a purely dielectric case for a reasonable dielectric contrast^. 



V. CONCLUSIONS 

Diamond and zinc-blende photonic crystals have been studied both in the purely dielec- 
tric case and in the presence of small inclusions of a low absorbing metal. In the former case 
we have shown that earlier plane wave method calculations of the photonic band structure 
were not converged. Unlike the case of a simple lattice (one scatterer per lattice primitive 
cell)@, for the case of a diamond lattice of dielectric spheres, even when using the plane-wave 
method based MIT ab initio program^!, one has to take a much higher number of plane waves 
than expected to reach a convergence comparable with the photonic KKR method. To reach 
convergence of the photonic band structure of a diamond lattice of dielectric spheres within 
1% the number of plane waves N c has to exceed 32 768 (cf. Ref.u) and still an extrapolation 
N — > ooi3 has to be performedB 

It has been shown that small inclusions of a low absorbing metal can have a dramatic ef- 
fect on the photonic band structure of diamond and zinc-blende structures. Several complete 
photonic band gaps (CPBG's) can open in the spectrum, between the 2nd-3rd, 5th-6th, and 
8th-9th bands. The respective metal-volume fractions to open a CPBG and to have a CPBG 
of 5% can be more than 40 times and 25 times smaller than in the case of an fee latticeSSEll. 
Unlike in the purely dielectric case, in the presence of small inclusions of a low absorbing 
metal and for a moderate dielectric constant (e < 10), the largest CPBG turns out to be the 
2nd-3rd CPBG. The 2nd-3rd CPBG is the most important CPBG, because it is the most 
stable against disorder!. One can create a sizeable (> 5%) and robust 2nd-3rd CPBG at the 
near infrared and visible wavelengths even for a photonic crystal which contains more than 
97% of low dielectric constant materials (e < 2.1). For diamond and zinc-blende structures 
of nonoverlapping dielectric and metallo-dielectric spheres, a CPBG begins to decrease with 
an increasing dielectric contrast roughly at the point where another CPBG starts to open 
- a kind of gap competition. A CPBG can even shrink to zero when the dielectric contrast 
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increases further (Figs. ||, |H], [EL]). As a consequence of the gap competition, the inclusions of 
a low absorbing metal have the biggest effect for the dielectric constant e G [2, 12], which is 
a typical dielectric constant at near infrared and in the visible for many materials, including 
semiconductors and polymers. Absorptance in the 2nd-3rd CPBG of 5% and located above 
A > 750 nm remains very small (< 2.6%) and is an order of magnitude smaller than in the 
case of an fee latticeEH Our zinc-blende structures require an almost twice as large f m as 
diamond structures do to achieve a comparable CPBG, yet they display smaller absorption 
due to an increased separation of the metal cores. The metallo-dielectric structures display 
a scalinglike behavior, which makes it possible to scale the CPBG from microwaves down to 
the ultraviolet wavelengths. Given the metal filling fraction f m , this scalinglike behavior is 
strongest for zinc-blende structures. 

Our results imply that any dielectric material can be used to fabricate a photonic crystal 
with a sizeable and robust CPBG in three dimensions, as long as metal inclusions can be 
added. These findings (i) open the door for any semiconductor and polymer material to be 
used as a genuine building block for the creation of photonic crystals with a CPBG and (ii) 
significantly increase the possibilities for experimentalists to realize a sizeable and robust 
CPBG at near infrared and in the visible. By relaxing the constraints on the dielectric 
constant, the ultimate goal of matching the photonic and electronic band gaps, which is a 
prerequisite for many applications!, can be achieved. 
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APPENDIX A: THE KKR METHOD 

Let A be a simple (Bravais) periodic lattice. According to the Bloch theorem, propagating 
wave tp in a periodic structure with the symmetry A is characterized by the Bloch momentum 
k. The latter describes translational properties of ip by any lattice vector r s G A, 

i/j(r + r s ) = ip(r) exp(zk • r s ). (Al) 

The Bloch property holds irrespective of the spin of a wave, i.e., is the same for scalar and 
vector waves. We shall confine ourselves to the case when, outside scatterers, wavefunction 
ip satisfies the scalar Helmholtz equation, 

[A + a 2 ]ip = 0, (A2) 

with a being a positive constant. Let Go\(a,k, R) denote the free-space periodic Green's 
function of the Helmholtz equation. The latter is defined as 

G 0A (a, k, R) = J2 G o(°, R " r s )e ik ' r * = ]T G (a, R + r^e^, (A3) 

r s GA r a eA 
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where R = r — r' and Go denotes a free-space scattering Green's function of the scalar 
Helmholtz equation at the points r and r'. In 3D 



G (a, R) = G (a, r, r') = _ e ^ wR \ (A4 ) 

Ann 

where R = |R|. Within the Korringa-Kohn-Rostocker (KKR) method^, band structure is 
determined by solving the KKR secular equation 

det [1 - t{a)g{cr,k)) = 0, (A5) 

where t is a single-scatterer T matrix and g is the matrix of structure constants^. Both 



t and g in Eq. (|A5| ) are considered as matrices with matrix elements labeled by pairs of 
angular momentum numbers (lm,l'm'), where —l<m<l. In the scalar case, for instance, 
in the case of multiple-scattering scattering of electrons, the matrix elements of g in the 
angular-momentum basis are defined as expansion coefficients of 

G 0A (cr,k,R) - Co((7,R) = ^W , rn'(o",k)j i ((xr)^ m (r)j / /((Tr / )F / r m ,(r / ), (A6) 

lm,l'm' 

where ji are the regular spherical Bessel functionsil, and Y\ m are the conventional spherical 
harmonicsii. 

It is interesting to note that, to a large extent, the scalar case also covers the scattering 
of vector and tensorial waves, i.e., waves with a nonzero soin, provided each field component 
independently obeys the scalar Helmholtz equation ( |A2| )Bii. In the vector case of multiple- 
scattering of electromagnetic waves, a = Uy/Fh/c, where u is the angular frequency, c is the 
speed of light in vacuum, and Eh is the dielectric constant of the host dielectric medium. Let 
L = (/m) stand for a multiindex of angular momentum numbers and let A label independent 
polarizations. The corresponding vector structure constants Gal,a' v are obtained from the 
scalar structure constants g^jj as 

1 i 

G A l,a>l> = E E U A (l,m,p,a)gi +P!m+a . j i, +pl}ml+al U AI (l',m l ,p',a'), (A7) 

p,p'=—l a,a'=—l 

where the Ua§ are group-theoretical coefficients, in the current case determined by the 
vector-coupling Clebsh-Gordon coefficientsffl. Similar for the acoustic and elastic waves. 
In the case of the layer KKR method§§. which can be viewed as an extension of dynamical 
x-ray diffraction theory^ to infinitely many diffracted beams, one is mainly interested in the 
reflection, transmission, and absorption properties of finite-width crystals. These properties 
are determined by the crystal scattering matrix. Roughly speaking, the denominator of the 
scattering matrix is the matrix which stands behind determinant sign in the KKR secular 
equation (|A5|), whereas the nominator of the scattering matrix is rather trivial. Therefore, 
in a numerical implementation, it is only required to make the following two modifications in 
the scalar KKR and LKKR numerical codes: (i) to include a single routine which performs 
the transformation of the scalar structure constants into vectorial structure constantsHlZH! 
and (ii) to use the routine which calculates t, the single-scatterer T matrix, appropriate to 
a given boundary value problem. The same is also true if one wants to adapt the scalar 
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3D KKR numerical codes which deal with clusters of more than 1000 particles of arbitrary 
shape and clusters of more than 1000 arbitrary impurities in a crystalSH. There seems 
to be no obstacle to have the photonic KKR methodsH§§ dealing with clusters of more 
than 300 particles. The KKR methods can be used for scatterers of arbitrary shapeEl and 
are optimized for lattices of spheres. 

If A is not a simple (Bravais) periodic lattice, i.e., there is more than one scatterer in the 
primitive lattice cell, the matrices t and g in Eq. ( |A5| ) become matrices with entries labeled 
by multiindices ALa, where A and L are as before and a runs over all the scatterers in the 
primitive lattice cell^j. 
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Figure captions 



Figure 1 - Photonic band structure for a close-packed diamond lattice of spheres of dielectric 
constant e = 12.96 in air (the same set of parameters as in RefJ). Two CPBG open in the 
spectrum with the respective gap to midgap frequency ratios of 1.3% and 4.2%. Frequency 
is in units [c/A], where c is the speed of light in vacuum and A is the lattice constant of 
a conventional cubic unit cell of the diamond lattice. If you prefer to use wavelength, the 
values on the y axis correspond to those of Aj A, where A is the vacuum wavelength. 
Figure 2 - Gap width to midgap frequency ratio for a diamond lattice of dielectric spheres as 
a function of the dielectric contrast e s /eh for sphere filling fractions f s = 0.17 and f s = 0.34. 
Figure 3 - Convergence of the band structure in Fig. [TJ as a function of the cutoff l max - 
Figure 4 - Photonic band structure for a close-packed zinc-blende lattice of spheres in air. 
One of the two spheres in the primitive lattice cell has a silver core and dielectric shell with 
refractive index n s = 1.45 (silica) with the core to total sphere radii ratio r c /r s = 0.75 
{fm — 7.2%), whereas the other is a homogeneous dielectric sphere with e = 12 of the same 
radius r s = 80 nm. Frequency is in units [c/A], as in Fig. [I]. 

Figure 5 - Calculated gap width to midgap frequency ratio of the 2nd-3rd CPBG for a 
close-packed diamond lattice of dielectric n s = 1.45 (silica) and n s = 2 (ZnS) coated silver 
spheres in air. The upper graph is for the sphere radius r s = 80 nm and the lower graph is 
for the sphere radius r s = 300 nm. The gap to midgap ratio is plotted as a function of the 
metal core radial filling fraction r c /r s . Metal volume fraction is then f m = 0.34 x (r c /r s ) 3 . 
Figure 6 - Calculated gap width to midgap frequency ratio of the 2nd-3rd CPBG for a close- 
packed metallo-dielectric zinc-blende lattice of spheres in air. Lattice primitive cell contains 
a silica sphere and a silver core - n s = 1.45 (silica) shell sphere of the same radii. g w is 
plotted as a function of the ratio r c /r s of the silver core - n s = 1.45 (silica) shell sphere for 
the cases r s = 80 nm and r s = 300 nm. Metal- volume fraction is then f m = 0.17 x (r c /r s ) 3 . 
Figure 7 - Calculated gap width to midgap frequency ratio of the 2nd-3rd CPBG for a 
zinc-blende lattice of spheres in air for the sphere radius of r s = 300 nm. One of the two 
spheres in the lattice primitive cell is a silver core - n s — 1.45 (silica) shell sphere, whereas 
the other is a dielectric sphere of the same radius. Figure shows the effect of the second 
sphere on the 2nd-3rd CPBG which is either a n c = 2 (ZnS) core - n s — 1.45 (silica) shell 
dielectric sphere with fixed r c /r s = 0.60 (dashed line) - or a homogeneous n = 1.45 (silica) 
sphere (solid line). 

Figure 8 - Gap width to midgap frequency ratio of the 2nd-3rd CPBG for a close-packed 
metallo-dielectric diamond lattice of silver core - dielectric shell spheres with r c /r s = 0.5 
{fm — 4.25%) in air as a function of the shell dielectric constant for the cases r s = 80 nm 
and r s = 300 nm. For a comparison, the dot-dashed line with diamonds shows the gap to 
midgap frequency ratio of the 2nd-3rd CPBG of a close-packed diamond lattice of purely 
dielectric spheres. 

Figure 9 - Gap width to midgap frequency ratio of the 2nd-3rd and 8th-9th CPBG's for a 
close-packed metallo-dielectric diamond lattice of silver core - dielectric shell spheres with 
r c/ r s — 0.5 (f m = 4.25%) in air as a function of the shell dielectric constant for r s = 300 nm. 
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For a comparison, the dot-dashed line with diamonds shows the gap to midgap frequency 
ratio of the 2nd-3rd CPBG of a close-packed diamond lattice of purely dielectric spheres. 
Figure 10 - Calculated gap width to midgap frequency ratio of the 2nd-3rd and the 5th-6th 
CPBG's for a close-packed metallo-dielectric zinc-blende lattice of spheres in air. One of the 
two spheres in the primitive lattice cell is a silver core - n s — 1.45 (silica) shell sphere with 
r c /r s = 0.75- whereas the other is a homogeneous dielectric sphere of the same radius r s 
(fm — 7.2%). The gap/midgap ratio is plotted as a function of the dielectric constant of the 
homogeneous sphere for the cases r s = 80 nm and r s = 300 nm. With increasing dielectric 
constant of the homogeneous sphere, first the 2nd-3rd CPBG opens. The appearance of the 
5th-6th CPBG roughly coincides with the point after which the 2nd-3rd CPBG begins to 
narrow. 

Figure 11 - Calculated gap width to midgap frequency ratio of the 2nd-3rd, 5th-6th, and 
8th-9th CPBG's for a close-packed metallo-dielectric zinc-blende lattice of spheres in air. 
One of the two spheres in the primitive lattice cell is a silver core - n s — 1.45 (silica) shell 
sphere with r c /r s = 0.75- whereas the other is a homogeneous dielectric sphere of the same 
radius r s = 80 nm (f m = 7.2%). The gap/midgap ratio is plotted as a function of the 
dielectric constant of the homogeneous sphere. Note that around e ~ 16 only the 5th-6th 
CPBG remains open. 

Figure 12 - Calculated gap width to midgap frequency ratio of the 2nd-3rd CPBG for 
a close-packed metallo-dielectric zinc-blende lattice of spheres in air as a function of the 
sphere radius. One of the two spheres in the primitive lattice cell is a silver core - n s — 1.45 
(silica) shell sphere with r c /r s = 0.75- whereas the other is a dielectric sphere of the same 
radius r s (f m = 7.2%). The dashed line corresponds to the case when the 2nd sphere is 
a homogeneous n = 2 (ZnS) sphere, the solid line is the case of a n c = 2 (ZnS) core - 
n s = 1.45 (silica) shell sphere with fixed r c /r s = 0.60, and the dotted line is for the case of 
a homogeneous n = 1.45 (silica) sphere. 

Figure 13 - An example of scaling of the midgap wavelength (solid line) of the 2nd-3rd 
CPBG for a close-packed zinc-blende lattice of spheres in air with the sphere radius. Gap 
edges are plotted by dashed lines. One of the two spheres in the lattice primitive cell is a 
silver core - n s — 1.45 (silica) shell sphere with fixed r c /r s = 0.75 - and the other sphere 
is a homogeneous n = 1.45 (silica) sphere of the same radius (f m = 7.2%). Gap width to 
midgap frequency ratio of the 2nd-3rd CPBG of the structure is shown by dotted line in 
Fig. 0. 

Figure 14 - An example of scaling of the midgap wavelength of the 2nd-3rd CPBG for 
a close-packed zinc-blende lattice of spheres in air with the sphere radius. One of the two 
spheres in the lattice primitive cell is n s = 1.45 (silica) shell sphere with a silver core 
with fixed r c /r s = 0.75, whereas the other is a purely dielectric sphere of the same radius 
(fm = 7.2%). The latter is either a homogeneous n = 1.45 (silica) sphere, a homogeneous 
n = 2 (ZnS) sphere, or a n c = 2 (ZnS) core-n s = 1.45 (silica) shell sphere with fixed 
r c/ r s = 0.6. Gap width to midgap frequency ratio of the 2nd-3rd CPBG of these structures 
is shown in Fig. (12] . 

Figure 15 - The extinction efficiency for a silver core - n s = 1.45 (silica) shell sphere with 
fixed r c /r s = 0.75 for different sphere radii. 

Figure 16 - An example of scaling of the midgap wavelength of the 2nd-3rd CPBG for a 
close-packed diamond lattice of coated silver spheres in air with the sphere radius. Spheres 
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are either n s = 1.45 (silica) coated silver spheres with fixed r c /r s = 0.6, or n s = 2 (ZnS) 
coated silver spheres with fixed r c /r s = 0.4. 

Figure 17 - Reflectance, transmittance, and absorptance of light incident normally on a 
two unit cells (12 planes) thick zinc-blende lattice of spheres in air stacked in the (111) 
direction. The spheres are as in Figs. [12] and [L3|. Dimensionless frequency is used on the x 
axis, where A is the lattice constant of a conventional unit cell of the cubic lattice. In all 
cases, the 2nd-3rd CPBG lies between « 0.7 and 0.8. Gap width to midgap frequency ratio 
of the 2nd-3rd CPBG of these structures is shown in Fig. the midgap wavelength of the 
2nd-3rd CPBG is shown in Figs. O and [TJ. 

Figure 18 - AABBCC stacking of our diamond and zinc-blende structures in the (111) 
direction. Stacking is along the body diagonal of the conventional unit cell of the cubic 
lattice. Lines perpendicular to the stacking z direction indicate hexagonal planes of spheres 
with the primitive lattice vectors &i and a 2 . Nearest-neighbor distance of hexagonal planes 
depicted by the same line is 1/3 of the body diagonal. Spheres in the planes separated 
by 1/4 of the body diagonal have identical xy positions. Spheres in the planes separated 
by 1/12 of the body diagonal are shifted by (ai + a 2 )/3 with respect to each other. The 
close-packed radius is 1/8 of the body diagonal. 

Figure 19 - Reflectance, transmittance, and absorptance of light incident normally on a two- 
unit-cell (12 planes) thick zinc-blende lattice of spheres in air stacked in the (111) direction 
as a function of frequency for geometrically identical cores made of different metals. The 
geometry of r s = 250 nm spheres is the same as in Fig. [17|. Dimensionless frequency is used 
on the x axis, where A is the lattice constant of a conventional unit cell of the cubic lattice. 
In all cases, the 2nd-3rd CPBG lies between 0.7 and 0.8. 



Figure 20 - The same as in Fig. [T^ except that the sphere radius is now 100 nm. 
Figure 21 - Behavior of the 2nd-3rd CPBG edges with an increasing relative radius of 
the silver core, r c /r s . Shown are the cases of a zinc-blende lattice having Ag@Si02 and 
Si02 spheres of the same radius r s = 80 nm in the primitive cell (solid lines), a diamond 
lattice of silver-core — n s = 1.45-shell spheres with radii r s = 80 nm (dot-dashed lines) and 
r s = 300 nm (dotted lines), and a diamond lattice of silver-core- n s = 2-shell spheres with 
radii r s = 80 nm (dashed lines) and r s = 300 nm (long-dashed lines). 

Figure 22 - Calculated gap width to midgap frequency ratio of the 2nd-3rd CPBG for a 
close-packed metallo-dielectric diamond lattice of spheres as a function of the host dielectric 
constant. The dependence is shown for the sphere radii r s = 80 and 300 nm. Spheres have 
a silver core and dielectric shell with e — 8. The core to total sphere radii ratio, r//, is 
either r c /r s = 0.5 (f m = 4.25%), or r c /r s = 0.6 (f m = 7.34%). For r c /r s = 0.4 (f m = 2.2%) 
and the dielectric shell with e = 4, the 2nd-3rd CPBG is closed in this range of the host 
dielectric constant. 

Figure 23 - Normalized axial force plotted against q/r s , the center offset from the focus in 
units of the sphere radius. Silver core —n s = 1.45 (silica) shell spheres with fixed r c /r s = 0.3 
are immersed in water = 1.33. Laser wavelength A increases by 100 nm from A = 600 
nm for r s = 80 nm to A = 1000 nm for r s = 200 nm. In all cases a stable equilibrium is 
achieved at around q/r s ~ 2. 
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